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Bipartite entanglement measures are fantastic tools to investigate quantum phases of correlated electrons. 
Here, I analyze the entanglement spectrum of gapped two-leg quantum Heisenberg ladders on a periodic ribbon 
partitioned into two identical periodic chains. Comparison of various entanglement entropies proposed in the 
literature is given. The entanglement spectrum is shown to closely reflect the low-energy gapless spectrum of 
each individual edge, for any sign of the exchange coupling constants. This extends the conjecture initially 
drawn for Fractional Quantum Hall systems to the field of quantum magnetism, stating a direct correspondence 
between the low-energy entanglement spectrum of a partitioned system and the true spectrum of the virtual 
edges. A mapping of the reduced density matrix to a thermodynamic density matrix is also proposed via the 
introduction of an effective temperature. 

PACS numbers: 75.10.Jm,05.30.-d,05.30.Rt 
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Introduction - The recent application of quantum informa- 
tion concepts to several domains of condensed matter [1] has 
proven to be extremely successful, giving new type of phys- 
ical insights on exotic quantum phases. Upon partitioning a 
many-body quantum system into two parts A and B, quantum 
entanglement can be characterized by the properties of the 
groundstate reduced density matrix of either one of the two 
parts, pa or ps- For example, entanglement entropies such 
as the Von Neumann entropy — Tt{pa ^Pa} or the family of 
Renyi entropies offer an extraordinary tool to identify a one- 
dimensional conformal invariant system [2] and provides e.g. 
a direct (numerical) calculation of its central charge [3]. 

Furthermore, the entanglement spectrum (ES) defined by 
the eigenvalues of a fictitious Hamiltonian H, where pa is 
written as exp {—%), has been shown to provide much more 
complete information on the system. In one dimension, un- 
derlying conformal field theory (CFT) leads to universal scal- 
ings of the ES (Ref. 4) and topological properties of the 
groundstate (GS) can be reflected by specific degeneracies [5]. 
Choosing a partition corresponding to a very non-local real- 
space cut, the ES has also been used to define non-local order 
in gapless spin chains [6] . 

Many-particle quantum entanglement is also a powerful 
tool to characterize topological features of two-dimensional 
GS (Ref. 7) as e.g. in dimer liquids on a cylinder geometry [8] . 
Also, bipartite ES have been shown to provide valuable infor- 
mations on the edge states of fractional quantum Hall states 
on spherical [9] and torus geometries [11] upon partition into 
two (identical) subsystems. Interestingly, the ES of the incom- 
pressible GS of a generic Landau-level-projected Coulomb 
Hamiltonian arranges into a low-energy CFT spectrum, a fin- 
gerprint of topological order, separated by an 'entanglement 
gap' from the high energy levels [9, 10]. 

Such advanced insightful analysis of the ES has not how- 
ever been fully exploited in low dimensional quantum mag- 
nets. In particular, the conjecture by Haldane of a precise cor- 
respondence between the entanglement spectrum and the true 
spectrum in reduced space, e.g. the spectrum of the subsystem 
A, is of very high interest and so far only supported by limited 
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FIG. 1: (Color online) (a) Ribbon made of two coupled periodic 
Heisenberg chains (2-leg ladder). The partition into two identical 
A and B subsystems is made by cutting the rungs along the dashed 
line, (b) Phase diagram of the two-leg ladder mapped onto a circle 
assuming J\ eg = cos 9 and J run g = sin 8. 



calculations on quantum Hall systems. [9, 11] Low dimen- 
sional quantum magnets offer a completely different class of 
many-body systems where new aspects of this correspondence 
can be investigated, giving further insights on this fascinating 
scenario. 

Model and System - In this manuscript, I consider a 2- 
leg ladder made of two quantum Heisenberg spin- 1/2 chains 
coupled via a "rung" exchange coupling J mn g, as shown in 
Fig. 1(a). Such a quantum magnetic ladder [12] offers an 
attractive although still simple system with three non-trivial 
phases, as shown in the phase diagram of Fig. 1(b), depending 
on the signs of the leg (i.e. within the chains) and rung Heisen- 
berg exchange couplings, parametrized as J\ es = cos 9 and 
•A-ung = sin 6 respectively. I shall not consider here the case 
when both couplings are ferromagnetic leading to a trivial 
fully polarized ferromagnet (lower-left quadrant). The physics 
of the other two phases (occupying the three remaining quad- 
rants) can be easily understood starting from the strong rung 
coupling limit, i.e. when J rung ~S> J\ eg . When J\ e& = spin 
singlets or triplets form on the rungs depending whether the 
rung coupling is antiferromagnetic (AFM) or ferromagnetic 
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FIG. 2: (Color online) Various S n entanglement entropies (n = 1 
Von Neumann Svn entropy, n = 2 and n = oo Renyi entropies) 
computed on 2 x L ladders of length up to L = 14, normalized by 
Lin 2 and plotted versus the angle 9. Note that £o = Soo = SW . 
Only the two right quadrants of Fig. 1(b) are considered. For com- 
parison, the fluctuation of Sz,A (normalized by L/4) is also plotted 
(stars). The corresponding ES are shown in Fig. 3 for the values of 6 
marked by arrows. Inset: effective inverse temperature jig (see text). 



(FM). For an AFM rang coupling J rung > 0, upon turning 
on a leg coupling of either sign the product of rang singlets 
smoothly evolves into a (gapped) 'rung singlet' phase. For a 
FM rung coupling ,/ run g < and a small AFM leg coupling 
Jieg > the ladder system can be mapped onto an effec- 
tive gapped spin-1 chain [24] yielding an effective 'Haldane 
phase' (Ref. 13, 14). Remarkably, such gapped phases remain 
stable all the way to the weakly coupled chain regime. The 
rung coupling is therefore a 'relevant' perturbation. For ex- 
ample, while the spectrum of the decoupled AFM chains sys- 
tem is the tensor product of two gapless Conformal Field The- 
ory (CFT) invariant low-energy spectra [25] of central charge 
c = 1, any finite J rung opens a gap. Note also that the two 
rung singlet phases for AFM and FM leg couplings labeled as 
(I) and (II) in Fig. 1(b) are smoothly connected to each other. 
Extension to frustrated inter-chain couplings is considered in 
EPAPS [23]. 

The finite size two-leg ladder of Fig. 1(a) is topologically 
equivalent to a ribbon which can be partitioned into two halves 
A and B preserving periodic boundary conditions. This of- 
fers a simple convenient setup to investigate the entanglement 
between the two chain subsystems as a function of their cou- 
pling J r ung- I report below the entanglement entropies as well 
as entanglement spectra in the two considered gapped phases, 
computed numerically on 2 x 10, 2 x 12 and 2 x 14 clusters. 
It is shown that the ES reflects the underlying CFT scaling be- 
havior of the isolated chains. This is remarkable, in particular 
in the strong rung coupling limit where the two subsystems are 
strongly entangled producing a short spin correlation length. 
Note that I am considering here a different setup than the one 



used by Kallin et al. [15] to calculate entanglement entropies 
on N-leg Heisenberg ladders. Indeed, in the latter case, the 
A subsystem was chosen to have a 'two-dimensional scaling' 
with its linear size. 

Results - Characterizing the entanglement between A and 
B requires the knowledge of the reduced density matrix pa 
of the A subsystem (i.e. the upper AFM chain). After com- 
puting the GS by Lanczos exact diagonalisation on finite 
2 x L periodic clusters, an explicit use of translation sym- 
metry enables to express pa in a block-diagonal form, where 
each block corresponds to an irreducible representation la- 
belled by one of the (allowed) total momentum K = 27r£, 
p = —L/2 + 1, ■ ■ ■ , L/2. These blocks can then be diago- 
nalised (separately) to compute the Von Neumann (VN) en- 
tropy, Svn = ~~ Tr{p.A In /3.a}, or the family of Renyi en- 
tropies, [16] S„ = j^lnTr{(pA) n }, n > 2. Note that 
SVn can De considered as lim n _>i S n = Si. Results for Svn 
and S2 in the Haldane and rung singlet phases (for Ji eg > 0) 
are reported in Fig. 2. The single-copy entanglement [17] ob- 
tained by taking the limit n — > 00 and given by S^ = — In Ao, 
where A is the largest eigenvalue of pA, is also shown for 
comparison. [26] An inspection of the finite size scaling of the 
data reveals that the leading term of all entanglement entropies 
is proportional to the size L (corresponding to the length of 
the edge between A and B) as expected from the area law. 
The data are therefore normalized by L In 2 which is the max- 
imum (exact) entanglement entropy obtained for the product 
of independent rung singlets (9 = n/2). The finite size cor- 
rections (details in EPAPS [23]) are found to be very small, 
almost not visible at this scale. As also expected, all S„ van- 
ish in the limit of decoupled chains, where the GS becomes 
a simple product state. Interestingly, the behaviors of Svn 
and S'2 are fairly similar, showing the same linear (quadratic) 
behavior with Ji eg ~ A9 in the strong rung coupling limit 
9 — > — 7r/2 (6> — > 7r/2). It should be noticed that, in contrast 
to Svn and S2, S x behaves linearly when 9 — > 7r/2, a behav- 
ior also seen in the quantum fluctuation of the (z-component 
of the) total spin Sz.a of the A subsystem. [18] 

I now move to the ES which is believed to contain much 
more information on the system. The ES is defined as the 
spectrum {S, a } of the hermician operator H given by the re- 
lation pa = exp {—%). The £ Q can then be obtained from 
the weights A Q of pa as £ Q = — In \ a . Typical ES (measured 
from the GS energy fo) plotted as a function of momentum 
K, for the three sizes, are shown in Fig. 3, both in the rung 
singlet (a,b) and the Haldane (c,d) phases. Note that the to- 
tal spin of the (A) subsystem is also a good quantum number 
which can be assigned to each level. It is remarkable that the 
low-energy excitations are spin-triplet that accurately resem- 
ble the des Cloiseaux-Pearson spectrum [19] of the quantum 
Heisenberg chain (up to a multiplicative factor); in particular 
two gapless modes at K = and K = n [27] are clearly 
visible. The lowest singlet excitations close to K = and 
K = n also form towers of states as predicted for the Heisen- 
berg chain [20]. This suggests strongly that the ES bears the 
same low-energy CFT structure. In that case, we expect, in 
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FIG. 3: (Color on line) Entanglement excitation spectra versus to- 
tal momenta K in the chain direction for four different values of 6 
(shown by arrows in Fig. 2) corresponding to the rung singlet (I) 
phase (a,b) and the Haldane phase (c,d). All low-energy excitations 
computed on 2 x 10, 2 x 12 and 2 x 14 ladders are shown by open 
(black) triangles, (red) squares and (blue) circles respectively. The 
lowest triplet eigenstates (for all L) are marked by (black) '+' sym- 
bols and are fitted as A£ = v\ sin (K)\ by dashed lines. The lowest 
singlet eigenstates for L = 14 are also marked by (red) ' x ' symbols. 



particular, the GS energy £ to scale as, 

Zo/L = e + d 1 /L 2 + 0(l/L 3 ). 



(1) 



Such a behavior is indeed found for strong AFM rung cou- 
plings where L 3> Z mag . Furthermore, the fit provides a num- 
ber in good agreement with the CFT prediction di = ncv/6, 
where v is the velocity of the triplet mode and the central 
charge is set to 1, e.g. d\ ~ 1.31 compared to ncv/6 ~ 1.24 
for 8 = jt/3. However, for smaller rung couplings at which 
L ~ /mag this scaling law is not satisfied (as expected) [23]. 

I finish the investigation of the entanglement spectra by 
considering the last case of the rung singlet (II) phase realized 
for a ferromagnetic leg coupling (Ji eg < 0) and an AFM rung 
coupling J rung > (upper left quadrant of the phase diagram 
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FIG. 4: (Color on line) (e,f) Same as Fig. 3 for the rung singlet (II) 
phase (only L = 10 and L = 14 are shown). Here the GS is the 
saturated ferromagnet. The total spin 5 of the lowest eigenstates are 
indicated by different symbols (and colors) and can be assigned to 
m-magnon bound states, m = 5 max — S = L/2 — S. The lowest 
energy excitations for L = 14 are fitted according to the formula 
for E m i n (K) (see text). In (e), the fit for L = 14 rescaled by a 
factor 14/10 (upper dotted line) gives also good agreement with the 
L = 10 data. 



of Fig. 1(b)) and smoothly connected to the limit of decoupled 
rung singlets (9 = n/2, Ji eg = 0). The results of the ES of 
2x10 and 2x14 ladders are shown in Fig. 4(e,f) for strong and 
weak rung couplings. At low energies, the ES are shown to 
coincide (up to an overall rescaling factor) with the spectrum 
of the ferromagnetic quantum Heisenberg chain, consisting of 
m-magnon bound states (or solitons) [22] given by E m (K) = 
2 J e ff sin 2 (K/2)/m, where J e g is an effective chain coupling. 
On a finite cluster, such multi-magnon excitations are subject 
to the kinematic constraint K > 2nm/L, where L is the lad- 
der length. Therefore, the lower-bound energy 'envelope' be- 
haves as, E min (K) ~ x^ eff s ^ n2 (K/2)/K up to small finite 
size corrections, as shown in Fig. 4 (note that E min (K) — > 
for all K in the thermodynamic limit). 

All these results on quantum ladders support the conjecture 
of a deep correspondence between the ES and the true spec- 
trum of the (virtual) edges. 

Effective temperature - Finally, I suggest that the reduced 
density matrix p& (in the GS) can be re-written as a thermo- 
dynamic density matrix by simply introducing an effective, 
model-parameter dependent, temperature scale Tg (focussing 
on the J^g > case). Indeed, comparison of Fig. 3(a) and 
Fig. 3(b) on one hand, and of Fig. 3(c) and Fig. 3(d) on the 
other hand, reveals almost identical spectra up to an overall 
multiplicative scale. This implies that pA can be written as, 
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where h is a parameter-free (extensive) Hamiltonian, 
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A 1 and /3e = T e 1 is an effective inverse temperature to be 
adjusted. Since, as shown in Fig. 3, the spectrum of h has 
the same c = 1 low-energy CFT structure as the Heisenberg 
chain Hamiltonian (up to a shift in GS energy) h can be 'nor- 
malized' by e.g. fixing the velocity v of the triplet branch 
to be UReis = t/2, the Heisenberg chain value. The effec- 
tive inverse temperature f3g is simply estimated from the ac- 
tual slope of the gapless (K = 0) mode of the correspond- 
ing ES and is reported in the inset of Fig. 2 as a function of 
6. Apart from logarithmic corrections, the thermal (magnetic) 
length is /id ~ Tg 1 (Ref. 21) which, heuristically, can be 
associated (up to a constant multiplicative factor of order 1) 
to the ladder correlation length Z mag . Therefore, the behav- 
ior of /3g in the inset of Fig. 2 simply reflects the behavior of 
Z mag with 6. In particular, in the strong AFM rung coupling 
regime J rung 3> J\ eg , f}$ is linear in 6 — n/2, in agreement 
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generally, within our normalization of h, /3g ~ / mag /2. Also, 
it is interesting here to use the exact equivalence between the 
entanglement VN entropy and the thermodynamic entropy of 
the (effective) finite-T subsystem. In the regime of weakly 
coupled chains, using the expression of the thermodynamic 
entropy of the Heisenberg AFM chain when T@ <C 1, one pre- 
dicts Svn/L ~ 7rTg/(3wHeis) (assuming again c = 1) which 
agrees (within less than 15% difference) with the calculated 
VN entropy, giving further support that h belongs to the same 
universality class as the AFM Heisenberg chain (for Ji eg > 0). 

Concluding remarks - In this paper, I showed that the ES 
of the (ground state) reduced density matrix of a 2-leg quan- 
tum ladder possesses remarkable universal features one can 
associate to its two single Heisenberg chain subsystems. This 
strongly supports a broader applicability (beyond quantum 
Hall systems) of the conjecture by Haldane establishing a deep 
correspondence between the ground state ES of a many-body 
system made of two entangled constituents with the true spec- 
tra of the virtual edges. For example, although for AFM leg 
coupling the two ground states of the quantum ladder at 9 > 
and 6 < belong to distinct topological sectors of the singlet 
spin Hilbert space [14] characterized by different 'string or- 
ders', it is remarkable that a similar c = 1 CFT low-energy ES 
is found and that Eq. (2) applies to both cases [28]. Similarly, 
although the two rung singlet phases (I) and (II) are smoothly 
connected, they exhibit completely different low-energy ES 
in straight connection to the different nature of their edges. 
Lastly, I notice that the results of this paper also apply to the 
case of frustrated ladders [23]. 
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This supplement of the main paper provides i) additional details and data regarding the scaling 
of the entanglement entropies and ii) extension of the calculations of ES to frustrated ladders. 

I. FINITE SIZE SCALINGS 

As reported in the main paper, finite size corrections of the entanglement entropies remain 
small. The finite size scaling of the Von Neumann entropy and of the single-copy entanglement 
are provided in Fig. 1 . In the regime of strong AFM rung coupling the magnetic correlation length 
'mag is very short and L » / mag . In this regime, the entanglement entropy scales with system size 

as, 

S VN /L = si + j- exp (-c 2 L// mag ) , (1) 

where the value of p is not crutial to the quality of the fit for L varying from 10 to 14 (typically 
p ~ 1 or ~ 2). However, finite size corrections at intermediate and small rung couplings do not 
show the exponential asymptotic behavior. Very likely, this is the sign of a cross-over regime for 
which L ~ / mag . Another possibility is that the subleading corrections in £ V n vs - L contain, in 
this regime, a constant term giving a 1/L contribution in Eq. (1). 

II. FRUSTRATED INTERCHAIN COUPLINGS 

Frustrated ladders as the one represented in Fig. 2(a) have been investigated in the literature for 
both antiferromagnetic interchain couplings J mng > and J x > [1], and ferromagnetic inter- 
chain couplings J r ung < and J x < [2]. A schematic phase diagram is shown in Fig. 2(b). In the 
two unfrustrated quadrants (not studied here) we recover the same rung singlet and Haldane phases 
as in the main paper. A cross-over is expected in the frustrated regions corresponding to the two 
other quadrants when 2 J x ~ J r un g - Weak coupling Renormalisation Group treatments [1, 2] pre- 
dict the occurrence of a very narrow strip [3] of a dimerized phase separated from the two limiting 
rung singlet and Haldane phases by continuous transitions (i.e. quantum critical points). Numer- 
ical evidences of this intermediate phase remain rather weak, except on the ferromagnetic side at 
large enough interchain couplings. Note also that, at large enough antiferromagnetic interchain 
couplings, a direct first-order transition has been found, as in the classical limit. 

The Von Neumann entropy as well as the n = oo Renyi entropy are shown in Fig. 3 in the 
two "frustrated" quadrants, for equal magnitudes of the exchange intra- and interchain couplings, 
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FIG. 1 : Finite size scaling of the Von Neumann entanglement entropy (left) and of the GS energy £o of the 
ES spectrum (identical to the single-copy entanglement Soo) (right), both normalized by system length L, 
as a function of 1/L 2 . Multiplicative factors have been used to bring all data in a single frame. The values 
of 9 are those indicated by the arrows in the Figure 2 of the main paper. The dashed ines correspond to an 
exponential fit in the left plot and to tentative power-law fits in the right plot. 



meaning one fixes ./J r 2 ung + J 2 , = J\ eg (i.e. = 7r/4). As in the main paper, the entanglement 
entropies have been normalized with respect to the ladder length L — 12 (and divided by In 2). 
The data clearly show a dramatic and rather sudden drop of the entanglement entropies in the 
intermediate regions 2 1 J x | ~ | J mng | . This is reminiscent of the behavior obtained in the main 
paper when the two chains get decoupled as J rung — > 0. The new remarkable result here is that 
the two chains effectively (almost) decouple although the magnitudes of the intra- and interchain 
couplings remain equal (i.e. */ J r 2 un „ + J 2 = J\ eg ) in the present example. On the other hand, the 



system size remains too small to reveal the presence of a finite intermediate region. However, this 
region, if it exists, should behave as very weakly coupled dimerized chains. Otherwise, it might 
correspond to a very exotic quantum critical point (for antiferromagnetic interchain couplings). 

Finally, I discuss briefly the ES of the frustrated ladders. The results shown here are fully con- 
sistent with the "conjecture" discussed at length in the main paper. In other words, the ES exhibit 
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FIG. 2: (a) Ribbon made of two coupled periodic Heisenberg antiferromagnetic chains (2-leg ladder) cou- 
pled via vertical J run g an d diagonal J x couplings. The partition into two identical A and B subsystems is 
made by cutting the rungs (and the diagonals) along the dashed line as in the main paper, (b) Phase diagram 
of the two-leg ladder for a fixed value of the AFM leg coupling Ji eg = cos 9 > 0. The parameter space 
( J m ng! 2 J x ) is mapped onto a circle of radius sin 9 i.e. J rung = sin 9 cos <j> and 2 J x = sin 9 sin <p. Frus- 
tration occurs when these interchain couplings are both FM or AFM. The thin shaded areas schematically 
represent the (narrow) maximally frustrated regions 2J X ~ Jrung where a dimerized phase is expected. 

a generic low-energy CFT structure (c = 1) similar to the one of an isolated AFM Heisenberg 
chain. The ES hence enables to characterize the "edges" of the frustrated ladder. 
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FIG. 3: Von Neumann entanglement entropy SVn entropy and single-copy entanglement 5^ (i.e. Renyi 
entropy for n = oo) computed on a 2 x 12 ladder, normalized by L In 2 and plotted versus the angle <f> (see 
Fig. 2(b)) for fixed 9 — n/4. Note that 5^ = £o ( see main paper). Only the two opposite quadrants of 
Fig. 2(b) exhibiting the narrow intermediate regions are shown. For convenience, S^ has been multiplied 
by a factor 2. Note the dip is more abrupt in the antiferromagnetic case. 



